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In this work we explore the background dynamics when dark energy is coupled to dark 
matter with a suitable interaction in the universe described by brane cosmology. Here 
DGP and the RSII brane models have been considered separately. Dark energy in the 
form of Generalized Cosmic Chaplygin gas is considered. A suitable interaction between 
dark energy and dark matter is considered in order to at least alleviate (if not solve) the 
cosmic coincidence problem. The dynamical system of equations is solved numerically and a 
stable scaling solution is obtained. A significant attempt towards the solution of the cosmic 
coincidence problem is taken. The statefindcr parameters arc also calculated to classify 
the dark energy models. Graphs and phase diagrams are drawn to study the variations 
of these parameters. It is also seen that the background dynamics of Generalized Cosmic 
Chaplygin gas is consistent with the late cosmic acceleration, but not without satisfying 
certain conditions. It has been shown that the universe in both the models follows the 
power law form of expansion around the critical point, which is consistent with the known 
results. Future singularities were studied and our models were declared totally free from any 
types of such singularities. Finally, some cosmographic parameters were also briefly studied. 
Our investigation led to the fact that although GCCG with a far lesser negative pressure 
compared to other DE models, can overcome the relatively weaker gravity of RS II brane, 
with the help of the negative brane tension, yet for the DGP brane model with much higher 
gravitation, the incompetency of GCCG is exposed, and it cannot produce the accelerating 
scenario until it reaches the phantom era. 
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INTRODUCTION 



Recent cosmological observations have indicated that the observable universe enters into an 
epoch of accelerated expansion . In the quest of finding a suitable model for universe, Cosmol- 
ogists started to investigate the root cause that is triggering this expansion. Within the framework 
of the general relativity, the acceleration can be phenomenally attributed to the existence of a 
mysterious negative pressure component which violates the strong energy condition i.e. p + 3p < 0. 
Because of its invisible nature this energy component is aptly termed as dark energy (DE) [3]. 
Moreover, quite surprisingly, observations spilled out definitively that about 70 percent of the Uni- 
verse is filled by this unknown ingredient and in addition that about 25 percent of this is composed 
by dark matter (DM). 

With the introduction of DE, search began for different candidates that can effectively play its 
role. DE represented by a scalar field 1 jj] is often called quintessence. Not only scalar field but 
also there are other Dark fluid models like Chaplygin gas which plays the role of DE very efficiently. 

get modified into Generalized Chaplygin gas 
In this context it is worth 



Extensive research saw Chaplygin gas (CG) [5j, 

u 



(GCG) |2HU| and then to Modified Chaplygin gas (MCG) 3, Q- 



mentioning that Interacting MCG in Loop quantum cosmology (LQCJwas studied by Jamil et al 
. Dynamics of MCG in Braneworld was studied by Rudra et al jl5|]. Other than these other 



forms of Chaplygin gas models have also been proposed such as Variab 



e Modified Chaplygin gas 

(VMCG) (16| and New Variable Modified Chaplygin gas (NVMCG) [l7]. Other existing forms of 
jis] ]. k-essence lpl. tachvonic field 



20JI, etc. 



DE are phantom 

In 2003, P. F. Gonzalez-Diaz [2l|] introduced the generalized cosmic Chaplygin gas (GCCG) 
model. The speciality of the model being that it can be made stable and free from unphysical 
behaviours even when the vacuum fluid satisfies the phantom energy condition. In the previous 
studies related to DE corresponding to phantom era Big-Rip is essential, as the time gradient of 
scale-factor blows to infinity in finite time. For the first time P. F. Gonzalez-Diaz including the 
GCCG model showed that Big Rip, i.e., singularity at a finite time is totally out of question. Hence 
in such models there is no requirement for evaporation of black hole to zero mass. The Equation 
of state (EoS) of the GCCG model is 



(1) 



1 in the presence of a scalar field the transition from a universe filled with matter to an exponentially expanding 
universe is justified 
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where C = — 1, with A being a constant that can take on both positive and negative values, 
and — C > 00 > 0, C being a positive definite constant, which can take on values larger than unity. 
GCCG can explain the evolution of the universe starting from the dust era to ACDM, radiation 
era, matter dominated quintessence and lastly phantom era 22|. In this context it should be stated 
that in 23f] Chowdhury and Rudra studied Interacting GCCG in Loop Quantum Cosmology (LQC). 

Currently, we live in a special epoch where the densities of DE and DM are comparable. Al- 
though they have evolved independently from different mass scales. This is known as the famous 
cosmic coincidence problem, till date several attempts have been made to find a solution to this 
problem 



24 



341 ] . A suitable interaction between DE and DM is required if we wish to find an 



effective solution to this problem. It is obvious that there has been a transition from a matter 
dominated universe to dark energy dominated universe, by exchange of energy at an appropriate 
rate. Now the expansion history of the universe as determined by the supernovae and CMB data 
31 , 32I bounds us to fix the decay rate such that it is proportional to the present day Hubble 



parameter. Keeping the fact in mind cosmologists all over the world have studied and proposed a 
variety of interacting DE models 3514411] . 



As we have stated earlier, modifying the right hand side of Einstein's equation (DE approach) 
was not the only way to explain the increase in the rate of the expansion. We can also modify the 
gravity part of the left hand side in order to demonstrate the present day universe. In this context 
Brane-gravity was introduced and brane cosmology was developed. A review on brane-gravity and 
its various applications with special attention to cosmology is available in 
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451 ] . In this work 



we consider the two most popular brane models, namely DGP and RS II branes. Our main aim 
of this work is to examine the nature of the different physical parameters for the universe around 
the stable critical points in two brane world models in presence of GCCG. Impact of any future 
singularity caused by the DE in brane world models will be studied. Finally some cosmographic 
parameters will be studied in brief. 

This paper is organized as follows: Section 2 comprises of the analysis in RS II brane model. 
Section 3 deals with the analysis in DGP brane model. In section 4, a detailed graphical analysis 
for the phase plane is done. In section 5, future singularities arising from our models are studied, 
followed by the study of some cosmographic parameters in section 6. Finally the paper ends with 
some concluding remarks in section 7. 
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II. MODEL 1: RS II BRANE MODEL 



Randall and Sundrum [4q, |47j] proposed a bulk-brane model to explain the higher dimensional 
theory, popularly known as RS II brane model. According to this model we live in a four di- 
mensional world (called 3-brane, a domain wall) which is embedded in a 5D space time (bulk). 
All matter fields are confined in the brane whereas gravity can only propagate in the bulk. The 
consistency of this brane model with the expanding universe has given popularity to this model of 
late in the field of cosmology. 

In RS II model the effectiv e equat ions of motion on the 3-brane embedded in 5D bulk having 
^-symmetry are given by 



43, 



47 



^Gpv = -A 4 q^ + kIt^u + rcflL^, - (2) 

where 



«4 = g Xl 4 > ( 3 ) 



A 4 = i f A 5 + | «|A 2 ) (J! 



2 ° V 6 
and 

!V = ~\ V* r ° + ^ TT ^ + g Q^T a ^T aP - — q^T 2 (5) 

and E^y is the electric part of the 5D Weyl tensor. Here k$, A5, t^ u and A4 are respectively the 
5D gravitational coupling constant, 5D cosmological constant, the brane tension (vacuum energy), 
brane energy- momentum tensor and effective 4D cosmological constant. The explicit form of the 
above modified Einstein equations in flat universe are 



K 

- n- + - 

and 



3H 2 = A 4 + 4 P+ -± p 2 + -^U (6) 



2H + 3H 2 = ^-^V-f x PV-f x P 2 -^U (7) 



The dark radiation U obeys 
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U + 4HU = (8) 

where p = p gC cg + Pm and p = p gC cg + Pm are the total energy density and pressure respectively. 

As in the present problem the interaction between DE and pressureless DM has been taken into 
account for interacting DE and DM the energy balance equation will be 

Pgccg + 3H (1 + UJgccg) Pgccg = ~Q, for GCCG and (9) 

p m + 3Hp m = Q, for the DM interacting with GCCG. (10) 

where Q = 3bHp is the interaction term, b is the coupling parameter (or transfer strength) and 
P = Pgccg + Pm is thetota 



p + 3H {p + p) = 
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cosmic energy density which satisfies the energy conservation equation 

152 1 . 

As we lack information about the fact, how does DE and DM interact so we are not able to 
estimate the interaction term from the first principles. However, the negativity of Q immediately 
implies the possibility of having negative DE in the early universe which is overruled by to the 
necessity of the second law of thermodynamics to be held 53]. Hence Q must be positive and 



small. From the observational data of 182 Gold type la supernova samples, CMB data from the 
three year WMAP survey and the baryonic acoustic oscillations from the Sloan Digital Sky Survey, 
it is estimated that the coupling parameter between DM and DE must be a small positive value 
(of the order of unity) , which satisfies the requirement for solving the cosmic coincidence problem 



and the second law of thermodynamics [54f | . Due to the underlying interaction, the beginning of 
the accelerated expansion is shifted to higher redshifts. The continuity equations for dark energy 
and dark matter are given in equations (9) and (10). Now we shall study the dynamical system 
assuming A4 = U = (in absence of cosmological constant and dark radiation). 

A. DYNAMICAL SYSTEM ANALYSIS 

Due to the complexity of the equations, it is very difficult to find direct solutions for this system. 
So in order to avoid these complex calculations, we undertake the dynamical system analysis for 
our further evaluations. In this subsection we plan to analyze the dynamical system. Before 
proceeding, the physical parameters are converted into some dimensionless form, given by 

x = \na, U = J§2> V = W 1 ^ 
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where the present value of the scale factor ao = 1 is assumed. Now using equations (1), (6), (7), 
(9), (10) and (11) we get, 



du „. , x do,, „ o / \ 9kIu (u + v) 

- = -36 <„ + „) + 3* - JtaoS" - 64- (« + *) + 2A(1 4 _1, (J,,) 



9k| (u + v) ; 



4A 2 (1-k| (u + v))' 



'AX 2 (1-k 2 4 (u + v)) 
k\ (u + v) 4 



2 2 



-C 



(12) 



d« , „ 2 , , , 3«fo (u + ^) 3 f / 4A 2 (l - «j (u + t;)) 2 n 2 N 

— = 36 (u + f )+3f — 3k 4 v (w + tu 7y— < C + 2 C 

dx y ' 4A 2 (l - k| (u + v)) u \ V k2(u + i>) , 



(13) 



Where, u^ccg 1 * s the EoS parameter for GCCG in RS II brane determined as 



RSII _ Vgccg 



4 (u + v) 



c + 



4AV (1 -Kj(u + v)f - 2Cnj (u + vf 



k\ (u + v) 4 



9CCg Pgccg 4A 2 U 2 (1 -kKu + v))' 

In the above calculations for mathematical simplicity we have considered a = 1. 



(14) 



1. CRITICAL POINTS 

The critical points of the above system are obtained by putting 4^ = = But due to 
the complexity of these equations, it is not possible to find a solution in terms of the involved 
parameters. So we find a numerical solution for the above system, by putting the following values 
to the different parameters appearing in the system. We take, 



6=1.5, K4 = l, A = 0.1, C = -l, 10 = -1 



and get the following critical point for the above dynamical system. 



u r = 1.51586 



v r = 1.79374 



(15) 
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The critical point correspond to the era dominated by DM and GCCG type DE. For the critical 
point (u c ,v c ), the equation of state parameter given by equation (14) of the interacting DE takes 
the form 



id 



RSII _ Pgccg 



K% (u c + V c f 



gccg 



Pgccg 4X 2 U 2 (l - k\ (u c + V c )) ' 



C + 



4A 2 u 2 (1 - k\ (u c + v c )) 2 - 2Ck\ (u c + v c 



K% (U c + V c ) 



(16) 



2. STABILITY AROUND CRITICAL POINT 

Now we check the stability of the dynamical system (eqs. (12) and (13)) about the critical 
point. In order to do this, we linearize the governing equations about the critical point i.e., 



u = u c + Su and v = v c + Sv, 



(17) 



Now if we assume / = ^ and g = ^ , then we may obtain 

6 (£) =l d uf}c^+[d v f] c Sv 



and 



' ^) = l d ^9] c Su + [d v g] c Sv 



(18) 



(19) 



where 



d u f = 



±k\u 2 (u + v) 2 (-1 + k\ {u + v)Y X 2 



Ck\ (u + vf {u (-3 + k\u) + v- k\v 2 } 



+2u 2 {-1 + k\ (u + v)} A |4/c| (« + u ) 4 (2u + u) A + 6u (it + 2u) A - 2k\ (u + v)(u(2-b + 6u) 
+ (2 - 6 + 12u) u) A + (u + w) 3 (3 (« + vf {3u + u) + 2 ((& - 10) u + {b - 6) v) a) + ^ (u + v) J 



(-3 (u + vf (4u + u) + 2 (« (8 - 26 + 3u) - 2 (b - 3 (1 + u)) v) a) } 



(20) 



d v f = - 



2k\u (u + v f (-1 + k\ {u + «)) 3 A 2 



C7e!| (u + vf {-2 + k\{u + v)} + An {-1 + k\ (u + v)} 



8 



| -6u 2 A + Ak\u (it + vf\ + 2k\ (u + v) (6u 2 + b(u + v)) \ + 2k\ (u + vf (su (u + vf + ((6 - 4) u + &u) a) 



+k 4 (u + u) 2 (-9u (u + vf - 2 (u (3u + 26 - 2) + 26t;) a) } 



(21) 



d u g 



4b - Ak\v - 



&k\vw (u + v) 2 ju + K 2 UV + v (— 1 + k 2 v) } < 


4« 2 (-l + K 2 (u+u)) 2 A 2 ' 


r 


{-1 + K 2 (U + V)} 


■\ 


Ck\ (u + v f - 4u 2 (-1 + k\ (u + v)f A 2 } 



2/^ (u + v)UC+[-C + 



4m 2 (-1+k|(u+i))) 2 A 5 
K t(«+f) 4 



(i-K 2 (u + v)yx 2 



/ r , , I n , 4« 2 (-l+ K |(u+t))) 2 A 2X 

3k 4 v(u + v) [C+[-C+ 4(u+v)i 



u{-l + K 2 (u + v)) 2 \ 2 

4v ( U + ,) 3 \c + (-c + 1 4{ :i f )} 



2 A 2X ~ W 



(-i + K 2 { u + v)y\ 2 



(22) 



a 3 
d v g = - A 



4 + 46 - 4k|w - Ak\ {u + v) + 



Sn 2 uv W (u + vf {4(u + v)-2}\-C + ^fgg^ ~ 



2/c^(u + ur < c+ -C + 



{k 2 (« + w ) - 1} \Ck\ (u + v) 4 - 4u 2 (k 2 (u + v)- l) 2 A 2 } 

4« 2 ( K 2 («+ 1 ;)-l) 2 A 2 



u 



{1-4 ( u +v)y\i 



3^(u + v) 2 lC+[-C + 



4M 2 (K 2 (n+«)-l) A 2 



2 , o \ ~ w 



U 



{-1 + K 2 (U + V)YX 2 



4 (u + vf { C + -C + 



4m 2 (k|(u+u) — l) A 2 
4(«+^) 4 



2,,\ —10 



7/ 



{-1 + K 2 (« + ?;)p 2 



(23) 
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The Jacobian matrix of the above system is given by, 



r(RSII) = 

\ Su 




The eigen values of the above matrix are calculated at the critical point (u c , v c ) and are found 
to be Ai = —1169.45, A2 = —3.38251. Hence it is a stable node. 

Description of figures 

Figs.l, 2, 3, 4 : The dimensionless density parameters are plotted against e- folding time for 
different values of brane tensions, A. The initial condition is v(0) = 0.9, u(0) = 0.2. The other 
parameters are fixed at w = —1, C = —1, K4 = 1 and b = 0.01. The brane tensions are respectively 
A = 1,-10, -50 and -100. 

Fig. 5 : The dimensionless density parameters are plotted against e- folding time for a high 
value of interaction (b = 5). The initial condition is v(0) = 0.9, u(0) = 0.2. The other parameters 
are fixed at w = — 1,C = —1, K4 = 1 and A = —10. 

Fig. 6: The dimensionless density parameters are plotted against e-folding time for same initial 
condition. The initial condition is t>(0) = 0.6,it(0) = 0.6. The other parameters are fixed at 
w = — 1,C = —1, K4 = 1 and A = —5. 

Figs. 7, 8 : The phase diagram of the parameters depicting an attractor solution are obtained 
for different values of brane tension. The initial conditions chosen are v(0) = 0.5, -u(0) = 0.6 
(green); v(0) = 0.6, it(0) = 0.6 (blue); v(0) = 0.7, n(0) = 0.6 (red); u(0) = 0.8, u(0) = 0.6 (brown). 
Other parameters are fixed at The other parameters are fixed at w = — 1,C = —1, K4 = 1 and 
b = 0.01. The brane tensions are respectively A = 1, —10. 

Fig. 9: The phase diagram of the parameters depicting an attractor solution are obtained for a 
high value of interaction (b = 1). The initial conditions chosen are v(0) = 0.5, u(0) = 0.6 (green); 
v(0) = 0.6, n(0) = 0.6 (blue); v(0) = 0.7, u(0) = 0.6 (red); v(0) = 0.8, u(0) = 0.6 (brown). Other 
parameters are fixed at The other parameters are fixed at w = —1, C = —1, K4 = 1 and A = — 10. 



Fig. 10 : The ratio of density parameters is shown against e-folding time. The initial conditions 
chosen are v(0)=0.6, u(0)=0.6. The other parameters are fixed at The other parameters are fixed 




Fig. 10 



Fig. 11 



Fig. 12 



Fig. 13 
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at w = —1, C = —1, A = 1, K4 = 1 and b = 0.01. 

Fig. 11 :The deceleration parameter is plotted against the EoS parameter. Other parameters 
are fixed at ip(RSii) = — 0.5. 

Figs. 12, 13 : The statefinder parameter r and s are plotted against the EoS parameter. Other 
parameters are fixed at ip(Rsn) = —0.5. 



3. NATURE OF COSMOLOGICAL PARAMETERS 
1. Deceleration Parameter: 

The deceleration parameter, q = — 1 — in this model is calculated as, 

f p ( (RSI I) Paccg 9 \ (-. (RSI I) Paccg \ \ 

" " + 2 (l + £) (24) 

We consider a dimensionless density parameter Q gccg = • in terms of this density parameter 
the expression for the deceleration parameter, q( RSII *> can be rewritten as, 



q (RSH) = _ 1+ i 7V—T\ ~ ( 25 ) 

It is evident that for A — > oo, we retrieve the result for the Einstein gravity, 



9£G = -1 + | (l + J™ g II] n gccg ) (26) 
Now since ft 9CC9 = = ^ and assuming ^ = V(hsjj) w e get, 

q{RSII) = _ i + 3 {(i + ^(rsii ) )4TJ i) ^ v + (i + 2^ rsii) )} 

2 (l+^(flS/J)) 
Now at the critical point, (u, i>) — >■ (n c ,v c ), and using equation (27) we get 



(RSn) a„ u „ {(^ + ^Rsn ) )4ccT ) ^r c +^ + ^(RSii))} 

q (RSH) = _ 1+ Z{RSII) , where Z (RSII) - ' 



2 (nall> (l + ^RSII)) 



(28) 
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Case I: 
If 



^(RSII) 



(Ugccg + 1) U + V 



(Ug CC g + 2)+2V 

then ^(ijs//) = 0. So we have q = —1, which confirms the accelerated expansion of the universe. 

Case II: 
If 

4>(rsii) = -i, 

then Z( RSII ^ — > — oo. Therefore q — > — oo. Hence in this case we have the super accelerated 
expansion of the universe. 

Note: In both the above cases we see that ip(RSii) < 0- We know that i^^rsii) = 2\- Since the 
energy density, p is always positive, therefore A should be negative (A < 0). Hence in order to 
realize the recent cosmic acceleration the RS II brane model should possess a negative 



brane tension. It is known from 43| standard model fields are confined on the negative tension 
(or visible) brane rather than the positive tension ("hidden") brane of the RS II model. Hence our 
result is consistent with the basic idea of formulation of the RS II brane model. 

In this scenario, the Hubble parameter can be obtained as, 

H = (29) 

^^{RSII) 1 

where the integration constant has been ignored. Integration of equation (29) yields 



a(t) = a t 3Z (Rsii) (30) 

which gives the power law form of expansion of the universe. In order to have an accelerated 
expansion of universe we must have < Z^sn) < § • Using this range of Z^sn) i n the equation 
qRSH __]__)_ ^Z(rsii)i i- e -> ec L n - (28), we get the range of qi RSH ^ as — 1 < qi RSH ^ < 0. This is 
again consistent with an accelerated expansion of the universe. 

2. Statefinder Parameters 

As so many cosmological models have been developed, so for discrimination between these 
contenders, Sahni et al [6s|] proposed a new geometrical diagnostic named the statefinder pair 
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{r,s}, where r is generated from the scale factor a and its derivatives with respect to the cosmic 
time t upto the third order and s is a simple combination of r and the deceleration parame ter q. 
Clear differences for the evolutionary trajectories in the r — s plane have been found 
statefmder parameters are defined as follows, 

a r — 1 



64 



fifj]. The 



(31) 



~ am ' 3(q - 1/2) ' 

The expressions for the statefmder pair (eqn.(31)) in the RS II model can be obtained in the form 



r { Rsii) = ( 1 " J (1 " 3Z (RSII) ) . (32) 

and 

S (RSII) = Z (RSII) (33) 

The trajectories in the r — s plane for various existing models can exhibit quite different be- 
haviours. The deviation of these trajectories from the (0, 1) point defines the distance of a given 
model from the ACDM model. The statefmder pair {r, s} can successfully differentiate between 
a wide variety of cosmological models including a cosmological constant, quintessence, Chaplygin 
gas, and interacting dark energy models. In a given model the pair {r, s} can be computed and 
the trajectory in the r — s plane can be drawn. Furthermore, the values of r, s can be extracted 



from future observations 



67, 



68]. 



Note: It is quite interesting to note that the pair {rmsil)i s (RSii)} yields the ACDM 

model {r E c s E g} = {1,0} when Z( RSII ) = 0. 



III. MODEL 2: DGP BRANE MODEL 



A simp 
model 



e and effective model of brane-gravity is the Dvali-Gabadadze-Porrati (DGP) braneworld 



55H57I] which models our 4-dimensional world as a FRW brane embedded in a 5-dimensional 
Minkowski bulk. It explains the origin of DE as the gravity on the brane leaking to the bulk at large 
scale. On the 4-dimensional brane the action of gravity is proportional to whereas in the bulk 

it is proportional to the corresponding quantity in 5-dimensions. The model is then characterized 

M 2 

by a cross over length scale r c = such that gravity is 4-dimensional theory at scales a « r c 
where matter behaves as pressureless dust, but gravity leaks out into the bulk at scales a >> r c 
and matter approaches the behaviour of a cosmological constant. Moreover it has been shown that 
the standard Friedmann cosmology can be firmly embedded in DGP brane. 
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It may be noted that in literature, standard DGP model has been generalized to (i) LDGP 
mode l by adding a cosmological constant 3], (ii) QDGP model by adding a quintessence per'ect 
fluid [59], (iii) CDGP model by Chaplygin gas [60] and (iv) SDGP by a scalar field [61]. In [62] 
the DGP model has been analysed by adding Holographic DE (HDE). 

While flat, homogeneous and isotropic brane is being considered, the Friedmann equation in 



DGP brane model 



5514571] is modified to the equation 



H 2 



P + — + — 

3 4r2 2r r 



(34) 



is 



M 

where H = ~ is the Hubble parameter, p is the total cosmic fluid energy density and r c = 
the cross-over scale which determines the transition from 4D to 5D behaviour and e = ±1 (choosing 
M 2 = 8itG = 1). For e = +1, we have standard DGP(+) model which is self accelerating model 
without any form of DE, and effective w is always non-phantom. However for e = —1, we have 
DGP(— ) model which does not self accelerate but requires DE on the brane. It experiences 5D 
gravitational modifications to its dynamics which effectively screen DE. Brane world scenario is 
actually a modified gravity theory. If we write the Einstein equation for brane world in terms of 
Einstein gravity then the extra term can be treated as the effective DE. But that is not the physical 
DE. Moreover this DE is applicable only in Einstein gravity. But here we will consider the physical 
DE in brane world. So we have introduced the GCCG type fluid in brane. 

Consequently using the Friedmann equation (34) and the conservation equations, we obtain the 
modified Raychaudhuri equation 



2H ) H 



-H{p + p) 



(35) 



A. DYNAMICAL SYSTEM ANALYSIS 

Just like the previous model, here also we proceed to perform the dynamical system analysis 
for DGP brane model: 

The system is obtained by using the equations (1), (9), (10), (11), (34) and(35) as given below, 



du „, . v „ /- nnp\ 2r\J (1 — u — v) 



9u (u + v ) 
rf. (1 — u — v) 1 



C + 



9u 2 



r\ (1 — u — v) £ 



C 



(36) 
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dv , . , 2rfv (1 — u— v) 4 

— = 3b(u + v)-3v+ c V / 
ax 3u(l + u + w) 



rf:(l — u — v) | Vt'cCI - u ~ v ) 



Where, oj®^ is the EoS parameter for GCCG in DGP brane determined as, 



u dgp = Pgccg = r*(l-u-v) 

Pgccg 9li^ 



c + < 

(1 — u — v) 



(37) 



(38) 



For mathematical simplicity, here also we have considered a = 1. 

Note: It is to be noted that e has been considered as —1 (DGP(-) model). Hence 
it does not produce a self-accelerating scenario of the brane model, but needs GCCG 
as dark energy (interacting) in order to realize the accelerating scenario. Thus the 
introduction of interacting GCCG in DGP brane is properly justified. 

1. CRITICAL POINTS 

The critical points for the above system (eqns. (36) and (37)) are calculated by putting ^ = 
= Here also due to highly complicated forms of the equations, it is difficult to get an explicit 
solution in terms of all the parameters. So we find the following solution in terms of the interaction 
parameter as below. 

u lc = 1 (l + VT^j , v lc = 1 (l - VT^ib) (39) 



u 2c = g (l - Vl - 4ft) , v 2c = g (l + Vl - 46) (40) 



The other variables are taken as: 



w = -1, C = -1, r c = 1, 

It is obvious from the above values that the critical point exists only for b < \. The critical point 
correspond to the era dominated by DM and GCCG type DE. For the critical point (uj c , v j c ), i = 
1, 2, the equation of state parameter (eqn. 38) of the interacting DE takes the form: 



DGP 

gccg 



(1 - u ic - v ic ) 
9ul 



\4 



c + 



9u 2 



r* (1 - u ic ~ vi 



C 



where 2 = 1,2 



2. STABILITY AROUND CRITICAL POINT 



Now if we write / = ^ and g = , then we can obtain the following expressions 



d u f 5u+ d v f 



Sv 



and 



where 
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(41) 



(42) 



(43) 



d u f = 



3u 2 (u + v + l) 2 |-9u 2 + Cr* (u + v - l) 4 } I ^ (u + v - 1) 



-C + 



9u 2 



-C 2 ^ (« + «-!)'' 



{3u 3 + („-!)(„ + I) 2 + u 2 (v + 3) - u (v + 1) (v + 5)} |-C + + 



9u 2 



Tc {u + V — 1) 



Cr 4 (u + v - 1) J 



jr 4 (u + v - l) 4 (3u 3 + («-!)(« + I) 2 + u 2 (v + 3) - u {v + 1) (u + 5)) + 9u 2 (-C + 



9u 2 



ri(u + v- 1) 



(u + u - 1) (it + v + l) 2 - 2 (2n 3 - u + v 3 + 2u 2 (u + 1) + u (v - 4) (v + 1))) } + 9u 2 {9u 2 (-C+ 



9u' 



r*(u + v-l) 



-j J (l - 2u - (u + v) 2 + b (u + u + l) 2 ) + r A c (u + v - if ((« - 1) (w + l) 2 (1 + 2w) 



+u 3 (3 + 2to) + u 2 (3 + v + 2w - 2vw) - u (v + 1) (5 + v + 2 (v + 1) w)) }] (44) 



*/- 5 



-96- 



{-' 



2 ^ ~ C + rUu+v-lf 



Y 



(u + v + l) 2 |-9u 2 + Cr^ (« + i; - 1) 4 | 



u 



j |c 2 rf (« + v - l) 7 (u + u 2 - uv - 2 (1 + v) 2 ) } 
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-C + 



9u< 



r*(u + v-l) 



+ Cr 4 (u + v- if jr 4 (u + v - l) 4 (u + u 2 - uv - 2 (v + l) 2 ) - 9u 2 (-3 + 2u 



W - (« + 3, „ - ^ (-a + ^ TT? ) } - * { (-c + 



9u 2 



ri( u + v- iy 



+r 4 (u + u - l) 3 (u - uv - 2 (v + l) 2 (w + l)+ u 2 (1 + 2w)^j } 



(45) 



3u 2 (u + v + l) 2 |-9u 2 + Crj (u + v - l) 4 } I r*(u + v-l) 



-C + 



9u l 



-2C 2 r*v {u + v-l) 7 



(1 + 2u 2 -v 2 + u(v + 5)) (-C + - 



9u 2 



c 4 (« + V ~ lf 



Cr 



4 {u + u - l) 3 (-2r 4 v (« + u - l) 4 (1 + 2u 2 



-v 2 + u (v + 5)) + 9u 2 (^-C + 



9u 2 



ri{u + v- 1) 



(b (u + v - 1) (u + v + l) 2 + 2v (2u (u + 3) + v + 



f / 9u 2 \ w 

v 2 )))+ 9u 2 { -9u 2 ( -C + j j (2v + 6 (u + v + l) 2 ) + 2r 4 v (u + t> - l) 3 (2u 2 (w + 1) 

I V rj(u + v — l)J\ ' 



(v 2 - 1) (2w + 1) + u (5 + u + 4w)) }] 



(46) 



<9„<? = 



2r^ (u + v -l) A \ -C + 



-9 + 9b- 



9u(u+v) 



■c + 



9u z 



ri(u+v-iy 



u (u + v + ly 



(- + » - ') 3 {-g + - (-g + ggrg) "} 

■u (u + -u + 1) 

2r < ( „ + „ _ 1)4 |_ c + _ (_ c + ^-5^) -} 



u(u + V + 1) 



\Sr{ 1 • 3c • u ^3 • Iw 



(-c + 



9u 2 



—1—w 



(u + v-l)(u + v + l) 



(47) 
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The Jacobian matrix of the above system is given by, 

/ Sf 5f\ 
j(DGP) = I 5H ^ 

(«.«) I Sg Sg 

\ 5u 5v ) 

Here we notice a very interesting feature of the model. We see that u c + v c = 1 for this model. Since 
the denominators of the above partial derivatives contain the term u + v — 1, hence they become 
indeterminate at the critical point, leading to a highly unstable scenario (chaos). As a result, 
determination of eigen values at the critical point is not possible for this model. So we resort to 
an alternative technique for our evaluations. We will consider a very small neighbourhood of the 
critical point, thus avoiding the chaos, and then try to calculate the eigen values at any convenient 
point in the neighbourhood sufficiently close to the critical point. This is purely based on the 
assumption that a sufficiently close neighbouring point will retain most of the properties of the 
critical point except the indeterminate nature. Our evaluations led us to the following eigen values 
for the given system: 

Ai = -2.97091, A 2 = 0.734899. Hence it is a saddle point. 
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Description of figures: 

Figs. 14, 15, 16, 17, 18 : The dimensionless density parameters are plotted against e-folding 
time for different values of interactions, b. The initial condition is v(0) = 0.9, u(0) = 0.2. The 
other parameters are fixed at w = — 1,C = —1 and r c = 10. The interactions are respectively 
b = 0.01,0.1,0.5, 1 and 10. 

Fig. 19: The dimensionless density parameters are plotted against e-folding time for same 
initial condition. The initial condition is v(0) = 0.6, u(0) = 0.6. The other parameters are fixed at 
w = -1, C = -1, r c = 10 and b = 0.01. 

Figs. 20, 21, 22 : The phase diagram of the parameters depicting an attractor solution are 
obtained for different values of interactions. The initial conditions chosen are v(0) = 0.5, u(0) = 0.6 
(green); v(0) = 0.6, u(0) = 0.6 (blue); u(0) = 0.7, u(0) = 0.6 (red); v(0) = 0.8, u(0) = 0.6 (brown). 
Other parameters are fixed at The other parameters are fixed at w = — 1,C = —1, r c = 10 and 
b = 0.01. The interactions are respectively b = 0.01, 0.1 and 1. 

Fig. 23 : The ratio of density parameters is shown against e-folding time. The initial conditions 
chosen are v(0)=0.6, u(0)=0.6. The other parameters are fixed at The other parameters are fixed 
at w = -1, C = -1, r c = 10 and b = 0.01. 

Fig. 24 :The deceleration parameter is plotted against the EoS parameter. Other parameters 
are fixed at ip(nsii) = —0.5. 

Figs. 25, 26 : The statefinder parameter r and s are plotted against the EoS parameter. Other 
parameters are fixed at b = 0.1, r c = 10 and a = 0.001. 

3. NATURE OF COSMOLOGICAL PARAMETERS 
1. Deceleration Parameter: 
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We calculate the deceleration parameter q = —1 — (H/H 2 ), in this model as, 

4r c 2 (l + 4cc,/ 



JDGP) 



y^fTSa ( y/Arl + 3a + e\/3a 



(48) 



where a = i . The above deceleration parameter can be written in terms of dimensionless density 
parameter r2 mcs = ^ ££ as 



JDGP) 



4r 2 f 1 + 4f C g P) 9Ccg 



V^r 2 + 3a y^r 2 + 3d + e\/3a 



(49) 



Now like the previous case we obtain $l mcg = ^^ 3 - = So from the previous equation we get, 



JDGP) 



-1 + 



4r 2 h+ w (DGP) ^ 



(50) 



V4r c 2 + 3a (y4r 2 + 3d + eV3aj 

We consider the first stable critical point. At the critical point (u,v) — > (u c ,v c ). Hence using 
equation (50) we get 



7c P) = -1 + ^Z {DGP) , w/iere Z (DGP) 



(51) 



2 { y/4r* + 3d ( V4r 2 + 3a + e V^) } ' 

Moreover we see that for a = 0, i.e., for p — > oo, we retrieve the results for Einstein gravity, as 



given below, 



(52) 



Considering the DGP(-) model we get, 



q (DGP) = _ 1+ z {DGP) , where Z {DGP) 



ArUl + U DGP) " c 



J 9ccg Uc+Vc 



(53) 



{ V4r 2 + 3a [y/Ar 2 c + 3a - v 7 ^) } ' 
If we do not want to spoil the successes of the ordinary cosmology, we have to assume the r c is 



of the order of the present Hubble scale Hq 1 [571 ]. Hence r c ^ 0. In the previous model we have 
seen that for Z(p GP ^ =0, q = —1 and the recent cosmic acceleration is effectively realized. From 
equation (53), we see that ^1 + U)g^ F ^ U ™+ Vc J = 0, since r c ^ 0. From this relation we get, 



(DGP) 
gccg 



U c + V c 
U c 



(54) 



Now since 



u c +v c 



> 1, we should have, 



(DGP) 



u gccg' ' < _1 



(55) 
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The above range for the EoS parameter indicates the phantom era for GCCG type 
DE. 

Hence it is evident that when Wg^ P ^ = — Uc ^ c , Z{dgp) = 0) an d hence q = — 1. This is 
consistent with the recent cosmic acceleration. 

Moreover the Hubble parameter can be obtained as, 

H = 2 , (56) 
where we have ignored the integration constant. Integration of eqn.(56) yields 



a(t) = a t 3 < DGP > , (57) 

which gives a power law form of the expansion. In order to realize the accelerating scenario of the 
universe, we should have zz^ DGP) > ^ *- e -' ^ < ^(dgp) < §■ Using this range of ^(dgp) i n the 
equation qi DGP ' > = — 1 + %Zm@py We get the range of qi DGP ^ as — 1 < qi DGP ^ < 0. Therefore 
the deceleration parameter is negative and hence the result is consistent with the fact that the 
universe is undergoing an accelerated expansion. 



2. Statefinder Parameters: 

In the DGP brane model, we have the following expressions for the statefinder parameters, r 
and s as given below, 

r (D GP) = (l " 3Z{D 2 GP) ^j (1 - 3Z {DGP) ) . (58) 

and 

S (DGP) = Z {DGP)- (59) 



IV. DETAILED GRAPHICAL STUDY OF PHASE PLANE ANALYSIS 

Figs. 1 to 5 and 14 to 18, shows the plots of density parameters u and v, respectively for RS II 
and DGP brane model. We see that in case of RS II brane as the brane tension A decreases (Figs. 1 
to 4), more and more irregularity creep in, as far as the DE density parameter u is concerned. But 
all the four figures show an energy dominated universe, consistent with observational data. In fig. 
5, with a larger value of interaction, we get a matter dominated universe (unphysical situation), 
which really indicates that the interaction coupling parameter, b should be a small positive value. 
In case of DGP brane, we see that with the increase in interaction between energy and matter, 
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the density parameters become more and more comparable to each other, giving a solution to the 
cosmic coincidence problem . An identical result was obtained for GCCG in LQC in [23j. Figs. 6 
and 19, shows almost the same results for identical scenario. In figs. 7, 8 and 9, Phase diagrams 
for RS II brane have been obtained. The figs. 7 and 8 have been generated for different values of 
brane tension A. It is seen that as the tension decreases there is a greater tendency of the solution 
moving towards an attractor, thus giving a perfect attractor solution. In fig. 9 with an higher 
interaction, we get a far better attractor solution, but the direction of flow is reversed. This can 
be attributed to the fact that as the interaction grows in magnitude DE interferes more and more 
with DM, until the matter loses its dominance and its place is taken by the energy, thus giving a 
perfectly energy dominated scenario. This shift of power may be responsible for many unexplained 
phenomena of cosmology including the present one. In the figs. 20, 21 and 22, we have obtained 
the phase diagrams for the DGP brane model, with gradually increased values of interaction. Just 
like the previous model we see that with the increase in interaction, there is a greater tendency of 
the flow going towards a specific attractor point. In the figs. 10 and 23, plots for the ratio of the 
density parameters have been generated against the e-folding time, for the two models respectively. 
In the plots, it is evident that ^ decreases with time, thus exhibiting an energy dominated scenario. 

Figs. 11 and 24 shows the plot of deceleration parameter, q against the EoS parameter, for the 
two models respectively. In the fig. 11, q remain in the negative level throughout the quintessence 
era thus exhibiting the recent cosmic acceleration for the RS II model. The only condition being 
the negativity of the brane tension, A. But in the fig. 24, we see that q remains in the negative 
level for u) < —1, i.e. only in the phantom region. This not only shows that GCCG is a DE fluid 
with far lesser negative pressure compared to other dark energy models, but also it shows that the 
combination of GCCG in DGP brane model gives an inferior model of the universe compared to 
the other combinations, like MCG in DGP brane (refer to [r]]). The statefinder parameters have 
been plotted in the figures 12 and 13 for RS II brane and in figs. 25 and 26 for the DGP brane. In 
the figs. 12 and 13, we see that r decreases whereas s increases with the increase in EoS parameter. 
But in case of DGP brane, r initially decrease, and then increase after reaching a minimum value. 
s increases with the increases in EoS parameter just like the previous model. 

V. STUDY OF FUTURE SINGULARITY 

We speculate that any energy dominated model of the universe undergoing an accelerated 
expansion will result in a future singularity. The study of dynamics of an accelerating universe in 
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the presence of DE and DM is in fact incomplete without the study of these singularities, which 
are the ultimate fate of the universe. It is known that the universe dominated by phantom energy 
ends with a future singularity known as Big Rip 69(, due to the violation of dominant energy 



condition (DEC). But other than this there are other types of singularities as well. Nojiri et al 
4] studied the various types of singularities that can result from a phantom energy dominated 
universe. These possible singularities are characterized by the growth of energy and curvature at 
the time of occurrence of the singularity. It is found that near the singularity quantum effects 
becomes very dominant which may alleviate or even prevent these singularities. So it is extremely 
necessary to study these singularities and classify them accordingly so that we can search for 
methods to eliminate them. The appearance of all four types of future singularities in coupled 
fluid dark energy, F(R) theory, modified Gauss-Bonnet gravity and modified F(R) Horava-Lifshitz 
gravity was demonstrated in [181 ] . The universal procedure for resolving such singularities that 
may lead to bad phenomenological consequences was proposed. In Rudra et al [15(] it has been 
shown that in case of MCG in Brane-world, both Type I and Type II singularities are possible. In 
Chowdhury et al [23[] it was shown that in case of GCCG in LQC, the universe is absolutely free 
from any type of singularities. We proceed to study the singularities for the present case: 



A. TYPE I Singularity (Big Rip singularity) 

If p — > oo , \p\ — > oo when a — > oo and t —> t s . Then the singularity formed is said to be the 
Type I singularity. 

In the present case by considering the GCCG equation of state from equation (1) we find that 
there is no possibility for Typel singularity, i.e., Big Rip singularity, since a > 0. This is in absolute 
accordance with P. F. Gonzalez-Diaz who has successfully shown that by considering GCCG as 
the DE, Big Rip can easily be avoided, thus giving a singularity free late universe. 



B. TYPE II Singularity (Sudden singularity) 

If p — > p s and p s ~ 0, then \p\ — > — oo for t — > t s and a — > a s , then the resulting singularity is 
called the Type II singularity. 

In this case we consider the equation of state for GCCG, like the previous case for our investi- 
gation. We see that if p — > p s and p s ~ 0, then \p\ — )■ for t — )■ t 3 and a — > a s . Hence there is no 
possibility of the type II singularity or the sudden singularity in case of GCCG, primarily because 
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a > and w < 0. 

C. TYPE III Singularity 

For t — > t s , a — > a s , p — > cc and \p\ — > oo. Then the resulting singularity is Type III singularity. 
It is quite evident from the equation of state of GCCG that it does not support this type of 
singularity. 

D. TYPE IV Singularity 

For t — > t s , a — > a s , p — > and \p\ — > 0. Then the resulting singularity is Type IV singularity. 
This type of singularity is not supported by GCCG type DE. 

As a remark, one should stress that our consideration is totally classical. Nevertheless, it is 
expected that quantum gravity effects may play significant role near the singularity. It is clear 
that such effects may contribute to the singularity occurrence or removal too. Unfortunately, due 
to the absence of a complete quantum gravity theory only preliminary estimations may be done. 

VI. CONSEQUENCES OF THE EXPANDING UNIVERSE: DISTANCE 
MEASUREMENT OF THE UNIVERSE 

Cosmography is a field where we are concerned with the measurement of the Universe. In fact 
there are many ways to specify the distance between two points. This is primarily because, in 
the expanding and accelerating Universe, the distances between co-moving objects are constantly 
changing, and Earth-bound observers look back in time as they look out in distance. The unifying 
aspect is that all distance measures somehow measure the separation between events on radial null 
trajectories, i.e., trajectories of photons which terminate at the observer. Here we will compute 
and discuss various cosmological distance measures such as the look-back time, luminosity distance, 
proper distance, angular diameter distance, co-moving volume, distance modulus and probability 
of intersecting objects. 

A. LOOK-BACK TIME 

As light travels with finite speed, it takes time for it to cover the distance related to the redshift 
it encountered. The difference between the age of the Universe now (at observation) and the age 
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of the Universe at the time the photons were emitted (according to the object) is defined as the 
Lookback time to an object. So, a look into space is always a look back in time. It is used to predict 
properties of high redshift objects with evolutionary models, such as passive stellar evolution for 
galaxies. Thus if a photon emitted by a source at the instant t and received at the time to then 
the photon travel time or the lookback time t$ — t is defined by [t3] 

f a da 

to~t= — 60 
Ja a 

where oo is the present value of the scale factor of the universe and can be obtained from (57) 
at t = tQ. The redshift is an important observable parameter as they can be measured easily from 
the spectral lines, and the redshift increases with the recession of the object from us. Look-back 
time is used to predict properties of high-redshift objects with evolutionary models, such as passive 
stellar evolution for galaxies. The redshift z can be defined by 



which gives the look-back time in the following form 



t-t = -^-\ (62) 



3ZH 



2 



For accelerating universe we have already get Z < | . Early universe is represented by z — > oo 
implies t — > and late universe z — > — 1, which equivalently implied t — > oo. Also z — > gives the 
present age t — > to of the universe. 



B. PROPER DISTANCE 

We know that light needs time to get from an object to the observer. Therefore we can define 
a distance that may be measured between the observer and the object with a ruler at the time 
the light was emitted, as the proper distance. When a photon emitted by a source at time tp and 



received by an observer at time t then the proper distance between them is defined by 



m 



; a ° da f to dt , . 

d = a — = a — (63) 
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which gives 

As far as the current epoch is concerned the proper distance may also called the co-moving 
distance (line of sight) of the Universe. So between two nearby objects in the Universe, the 
distance between them remains constant with epoch if the two objects are moving with the Hubble 
flow. In other words, it is the distance between them which would be measured with rulers at the 
time they are being observed divided by the ratio of the scale factor of the Universe (a at that time 
to a now). Next thing to be defined is the transverse co-moving distance, which is a quantity used 
to get the co- moving distance perpendicular to the line of sight. For flat universe, the transverse 
co- moving distance is always identical to the co- moving distance (line of sight). That means the 
transverse co-moving distance = proper distance = d, for any model following the power law form 
of expansion. 



C. LUMINOSITY DISTANCE 



If L be the total energy emitted by the source per unit time and I be the apparent luminosity 
of the object then the luminosity distance is defined by 



70, 



7J 



d(l + z) 



H (3Z - 2) 



(1 + *) 



3Z 

(1 + z)- 



(65) 



D. ANGULAR DIAMETER DISTANCE 



We define the angular diameter of a light source of proper distance D observed at to by 

s D(l + z) 2 



70, 



d L 



(66) 



Now the ratio of the source diameter to its angular diameter (in radians) is defined as the angular 
diameter distance dA as furnished below, 



d A = - = d L (l+zy 2 = d(l + z)- 1 



(67) 
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Fig. 27 



For our models following power law form of expansion angular diameter distance id a) is given 



by, 



d A 



1 



1 



(68) 



H (3Z-2) \ l + z + 

It is used to convert angular separations in telescope images into proper separations at the 
source. It is famous for not increasing indefinitely as z —> oo; it gets inverted at z ~ 1 and 
thereafter more distant objects actually appear larger in angular size. The angular diameter 
distance is maximum at 



3(2-Z) 



1 



3Z - 2, 

and corresponding maximum angular diameter d,A\max taking the form 



Aa\i 



1 



#0 (3Z - 2) 



2 1+ 3(2-Z) 



1 



3Z-4 



3(Z-2) 



2 I A_W=z) 



1 



3Z-4 



3(Z-2) 



(69) 



(70) 



All the above four parameters have been plotted against the redshift parameter, z in fig. 27. 
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E. CO-MOVING VOLUME 



The co-moving volume Vc is the volume measure in which number densities of non-evolving 



objects locked into Hubble flow are constant with redshift. We define it as 

dV c = D H ^ + Z ) dA dQdz = -L (l + zf-'^ d\dVLdz (71) 

where is the solid angle element and d& is the angular diameter, E{z) = -fflf and Dh = -g^ 
is the Hubble distance (c is the velocity of light) and in our model we assume c = l,d = 1, Hq = 
72km/s/Mpc. 

So, the co-moving volume is proper volume times the ratio of scale factors now to then to the 
third power. Co-moving volume element are drawn in figure 28. We see that there is a gradual 
increase with increase in redshift z. 



F. DISTANCE MODULUS 



The distance modulus is define by 



D " = 5l0g (^) (72) 

because it is the magnitude of difference between objects observed bolometric (i.e., integrated 
over all frequencies) flux and what it would be if it were at 10 pc (this was once thought to be 
the distance to Vega) and di is the luminosity distance. Distance modulus Dm as a function of 
redshift have been shown in figure 29. We see that for z > 0, we do not get any plot for Dm- But 
for z < 0, Dm increases as z decreases. 



G. PROBABILITY OF INTERSECTING OBJECTS 



It is defined as the incremental probability dP that a line of sight will intersect one of the 
objects in redshift interval dz at redshift z. It is given by 



70, 



/i _i_ _A 2 

dP = n(z)a(z)D H [ Z) dz 
E{z) 



(73) 



where n(z) is the co-moving number density and <j(z) areal cross-section. Assuming n{z)a{z) = 1, 
we obtain 



1 , > o 3Z 

dP = — (1 + z) 2 '— dz 



(74) 
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For our model the expression of Probability of intersecting objects becomes 

2 



{< 



1 + z) 



1 



(75) 



H (6 - 3Z) 

In figure 30 we draw intersection probability P as a function of redshift. We see that, P increases 
as z increases. 



VII. CONCLUSION 



In this work, we have considered a combination of Generalized cosmic Chaplygin gas in standard 
Brane-world models. Two different models, namely the RS II brane and the DGP brane models 
have been considered for our evaluations. Our basic idea was to study the background dynamics 
of GCCG in detail when it is incorporated in brane gravity. Because of the complexity of the 
expressions it was impossible to find direct solutions for the system. So we resorted to dynamical 
system analysis for our computations. Dynamical system analysis was successfully carried out, 
critical points were found and the stability of the system around those critical points was tested. 
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Graphical analysis was done to get an explicit picture of the outcome of the work. In order to 
find a solution for the cosmic coincidence problem, a suitable interaction between DE and DM 
was considered. Figures of density parameters were drawn for different values of interaction. It 
was found that increase in interaction resulted in more and more comparable values of the density 
parameters of GCCG and DM. Since the tendency of DE domination over DM is lesser in case 
of GCCG compared to MCG, GCCG is identified as a dark fluid with a lesser negative pressure 
compared to MCG or any other forms of DE. 

It was found that GCCG in RS II brane is consistent with the late cosmic acceleration only if 
the brane tension is negative. For GCCG in DGP brane an accelerated expansion is realized only 
in the phantom era of the DE. In the quintessence era there is no possibility of an accelerating 
scenario. This is a very important result as far as modern cosmology is concerned. This really 
shows the less effectiveness of GCCG as a DE compared to others, like MCG which produced an 
accelerating scenario in DGP brane in the quintessence era itself. From the above results we 
can come to the conclusion that although GCCG with a far lesser negative pressure 
compared to other DE models, can overcome the relatively weaker gravity of RS II 
brane, with the help of the negative brane tension, yet for the DGP brane model 
with much higher gravitation, the incompetency of GCCG is exposed, and it cannot 
produce the accelerating scenario until and unless it reaches the phantom era. 

The dynamical system of equations characterizing the system was formed and a stable scaling 
solution was obtained. Hence this work can be considered to be a significant one, as far as the 
solution of cosmic coincidence problem is concerned. Study of future singularities had been carried 
out in detail. The model was investigated for all possible types of future singularities. From 
the above analysis we conclude that the combination of GCCG in brane gravity gives a perfect 
singularity free model (just like GCCG in LQC) for an expanding universe undergoing a late 
acceleration. Statefinder parameters were calculated and plots were generated for them. The 
evolutionary trajectories obtained, when compared with those of the ACDM model gave clear 
differences, thus characterizing the GCCG type DE irrespective of the theory of gravity. The 
deceleration parameter was also calculated and plotted against the EoS parameter. From these plots 
the notion of an accelerating universe was clearly realized. Finally some cosmographic parameters, 
involving different types of distance measurements have been studied for both the models, following 
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the power law form of expansion and plot were generated to characterize the parameters. 
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